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Abstract. - The three-spined stickleback Gasterosteus spp. is a model organism largely used in biology. Four 
species have been described in Europe: G. aculeatus , G. gymnurus , G. islandicus and the extinct G. crenobion- 
tus. Our integrative taxonomy study, including 194 specimens, demonstrates that these taxa are distinguishable 
neither with the mitochondrial COI marker nor the morphological characters like lateral plates, supposed to dif¬ 
ferentiate them. Thus, we invalidate Gasterosteus gymnurus , and we consider it as junior synonym of Gasteros¬ 
teus aculeatus. The taxonomical status of Gasterosteus islandicus still needs to be clarified. 


Resume. - Invalidity de Gasterosteus gymnurus (Cuvier, 1829) (Actinopterygii, Gasterosteidae) confirmee par 
la taxonomie integrative. 
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integrative, comprenant 194 specimens, demontre que ces taxons ne peuvent etre distingues ni avec le marqueur 
mitochondrial COI, ni avec les donnees morphologiques comme les plaques laterales. Nous invalidons done 
l’espece Gasterosteus gymnurus , et nous la considerons comme synonyme junior de Gasterosteus aculeatus. Le 
statut taxonomique de Gasterosteus islandicus necessite d’etre clarifie. 


The three-spined stickleback 
Gasterosteus aculeatus Linnaeus, 
1758 (Actinopterygii, Gasterostei¬ 
dae) has for decades been a model 
organism for evolutionary, behav¬ 
ioural, developmental and ecological research, which quali¬ 
fied it as “supermodel” (Merila, 2013). It inhabits sea or 
freshwater, and has a very large distribution in circumarctic 
and temperate regions (Froese and Pauly, 2014). Its taxono¬ 
my has often been discussed as it is very polymorphic, espe¬ 
cially regarding plate morphs. 47 species were described 
throughout the world (Eschmeyer, 2014), and have been 
mainly considered to be synonyms to the species present in 
Western Europe, G. aculeatus and Gasterosteus gymnurus 
Cuvier, 1829 (Kottelat, 1997). Four other species are cur¬ 
rently considered valid: Gasterosteus wheatlandi Putnam, 
1867, Gasterosteus islandicus Sauvage, 1874, Gasteros¬ 
teus crenobiontus Bacescu and Mayer, 1956 (extinct), and 
Gasterosteus nipponicus Higuchi et al., 2014, respectively 
endemic to North-West America, Iceland, Central of Europe 
and Japan. Considering morphological and molecular studies 


( e.g. Miinzing, 1962; Makinen and Merila, 2008), Kottelat 
and Freyhof (2007) changed the assignation of West Euro¬ 
pean sticklebacks into G. gymnurus, because of the low plate 
numbers on the flanks and the absence of a keel on the cau¬ 
dal peduncle. Keith et al. (2011) considered it to be the only 
species of three-spined stickleback in France, whereas Igle- 
sias (2012) recognized both species. However, several stick¬ 
lebacks “forms” in France, especially in the coastal basins 
of the Channel (North-West of France), display a variabil¬ 
ity in their lateral plate numbers, as already noted by Bertin 
(1925). This is of particular importance, as G. gymnurus was 
described solely on the character “presence of plates only 
on pectoral region” (translated from French: “G. gymnurus 
(...) n’a de ces plaques que dans la region pectorale”; Cuvier 
1829). A morphological and morphometric study (Wolt- 
mann and Berg, 2013) on more than 5000 German individu¬ 
als belonging to both forms found no difference between 
G. aculeatus and G. gymnurus. 

The aim of this study is to verify the status of G. aculea¬ 
tus and G. gymnurus using DNA taxonomy s ensu Tautz et al. 
(2003), combining morphological characters and the mito- 
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chondrial gene coding for cytochrome C oxidase subunit 1 
(COI) on a large circumpolar sampling, in order to evalu¬ 
ate both species in an integrative taxonomy approach (Padial 
et al., 2010). It occurs in a context where the taxonomy of 
European ichthyofauna undergoes many discoveries and 
changes (Kottelat and Freyhof, 2007; Denys et al., 2014; 
Geiger etal., 2014). 

MATERIALS AND METHODS 

Sampling 

Samples were collected between 2003 and 2013 with col¬ 
laborations of the FREDIE program (http://www.fredie.eu), 
the French National Agency for Water and Aquatic Environ¬ 
ments (Onema), and the Federation nationale de la peche 
en France (FNPF). A total of 174 specimens were caught, 
mainly by electrofishing, in 30 locations in the major French 
drainages. Ninety-five individuals from 37 other European 
sites were added (Fig. 1; Annexe 1). Specimens were fixed 
in 95% EtOH for later fin clipping. We did all morphological 
identification (juveniles excepted): G. aculeatus has a com¬ 
plete series of 29-35 bony scutes covering the trunk and the 
caudal peduncle (Fig. 2A), whereas G. gymnurus has only 
2-10 scutes on the anterior part of the trunk (Fig. 2B), and 
G. islandicus has a deep notch on the anterior margin of the 
pelvic girdle (vs. straight or rounded for the other taxa) (Kot¬ 
telat and Freyhof, 2007). One hundred sixty two specimens 
were used for DNA analyses. 

DNA extraction, PCR, sequencing and quality con¬ 
trol were performed according to Geiger et al. (2014) and 
for French samples Denys et al. (2014), yielding a dataset 


of partial COI sequences (627 bps). It was completed by 
sequences from other Gasterosteus species with locality 
and photo vouchers available on the Barcode of Life data¬ 
base (BOLD, www.boldsystems.org; Ratnasingham and 
Hebert, 2007) and GenBank. Several other gasterosteid spe¬ 
cies were used as outgroups (Fig. 1; Annexe 1). All new COI 
sequences were deposited in the EUEPI project folder in 
the BOLD with their voucher information. Alignment was 
performed manually, as neither marker includes indels. Phy¬ 
logenetic analyses were performed with Bayesian inference 
(MrBayes 3.2, Ronquist et al., 2012), with the GTR+I+G 
model selected by JModelTest 2.1.1 (Darriba et al., 2012). 
Two runs of two analyses with 10 million generations and 
sampling every 200 generations were performed, and 10% 
of trees were eliminated as burnin after checking for conver¬ 
gence. Intra- and inter-specific distances (p-distances) were 
calculated with MEGA 6 (Tamura et al., 2013). 


RESULTS 

The morphological identification according to Kottelat 
and Freyhof (2007) of the 159 non-juvenile three-spined 
sticklebacks (i.e. over 30 mm SL) placed G. islandicus (one 
population, four specimens) in Iceland, G. gymnurus (40 
populations, 100 specimens) from Portugal to Turkey includ¬ 
ing United Kingdom and Mediterranean drainages, and 
G. aculeatus (28 populations, 51 specimens) from Channel 
coast and east of Europe to Alaska and Canada. Mediterra¬ 
nean and all French populations were determined as G. gym¬ 
nurus, except one from Taute stream in Normandy (MNHN 
2013-1299), which was identified as G. aculeatus. The pop- 


Figure 1. - Geographical distribution of the vouchers 
used in the analysis in polar view. Gasterosteus acu¬ 
leatus (black squares), G. gymnurus (white squares), 
G. islandicus (white circle), G. wheatlandi (white stars) 
and unidentified morphologically specimens (grey 
squares). Asterisks mean approximative locations. 



38 


Cybium 2015,39(1) 








Denys etal. 


Invalidity o/Gasterosteus gymnurus 


ulations from the North Sea (ZFMK:ICH:55427), Acheron 
in Greece (ZFMK:ICH:55875 to 55879), and the Japanese 
specimen (UW 41883) were also identified as G. aculeatus. 

Genetic analyses were performed on 194 individuals 
for the COI marker. Twenty-four specimens belonging to 
G. aculeatus have identical sequences to specimens identi¬ 
fied as G. gymnurus. 

In the Bayesian phylogenetic tree reconstruction (Fig. 2), 
a deep divergence (12.62% mean pairwise divergence) sepa¬ 
rates G. wheatlandi from the strongly supported (1 ppv) 
clade including G. aculeatus , G. gymnurus and G. islandi- 
cus. The clade is a large polytomy; the variability within it is 
0.36% pairwise divergence, and it has no taxonomical struc¬ 
ture. Two G. gymnurus individuals from the United King¬ 
dom (North Sea: Ex70G2 and Ex70G4) are included into a 
Canadian subclade of G. aculeatus supported by one synapo- 
morphy (A in position 401). However, several subclades are 
geographically homogeneous, like a Western Mediterranean 
subclade from France to Italy (except for Brague stream) 
supported by two synapomorphies (G in positions 14 and 
547), with moderate support (0.87 ppv) and mean divergence 
from the others of 0.5%. On the other hand, some localities 
from surrounding areas or even areas within the geographic 
areas of these subclades, have different haplotypes. 

DISCUSSION 

Miinzing (1962) drew a map of Europe with the reparti¬ 
tion of the different plate morphs. He distinguished clearly 
the plateless morph “ leiurus ” in the West part of Europe, 
whereas the morph with plates “ trachurus ” occurs from 
Channel to Baltic Sea. Our morphological identifications 
demonstrate globally the same pattern. But one population 
in Channel drainage is fully plated; other French and Eng¬ 
lish populations have just the lower part of the flank covered. 


Bertin (1925) noticed that “forms” covered by plates “ tra - 
chura ” occurred only in costal basin of Channel, but dur¬ 
ing their growth, they passed by the other forms having less 
plates: “hologymna ”, “gymnura” and “ semiarmata ” (see 
Bakker and Sevenster, 1988 for the terminology review). 
This is why care must be taken to compare only adults, which 
is the case of our specimens. Moreover, one Greek popula¬ 
tion (Acheron drainage) is fully plated, whereas all other 
Mediterranean populations have a few plates. Thus, the mor¬ 
phological identification demonstrates that the species dis¬ 
tinction between G. aculeatus and G. gymnurus using plate 
numbers and arrangement as key character does not lead to a 
clear biogeographic pattern. 

The lateral plate number is controlled by the ectodyspla- 
sin-A Eda gene, with a stronger expression in adults than 
juveniles (see Barrett, 2010). Plate presence is strongly cor¬ 
related to the habitat: marine specimens are totally covered 
with plates, whereas freshwater populations are less shielded 
(Barrett 2010; Merila 2013). A coastal-inland gradient as 
positive relationship between salinity and plate numbers 
was documented (Raeymaekers et al ., 2007), activating 
firstly the anterior plates development, followed by posterior 
plates and finally middle plates (Bell, 2001). So plate reduc¬ 
tion probably constitutes an adaptation to freshwater habitat 
(Le Rouzic et al., 2011) colonized by fully plated ancestral 
marine populations, and often occurring in numerous cases 
through allelic substitution in the Eda gene (Leinonen et al., 
2012). The main explanation for these phenotypes is that 
lateral plates serve as defensive structures against the teeth 
of predatory fish in marine environments, but they alter the 
escape capacity in freshwater (Barrett, 2010). 

Density of plates is therefore not an appropriate taxo¬ 
nomical criterion for European three-spined sticklebacks, 
but represents an intraspecific polymorphism induced by 
habitat, contrary to G. nipponicus , which has a characteristic 
plate morph and is supported by ecological, chromosomal 




Figure 2. - A: Gasterosteus aculeatus , 
MNHN 2013-1299, FFFtagl2285, 
35 mm SL, Taute (Douve drainage) 
at Tribehou, France, 17 Sep. 2013; 
B: G. gymnurus , MNHN 2014-0010, 
FFFtagl2314,42 mm SL, Blaise (Seine 
drainage) at Saint-Ange-et-Torgais, 
France, 24 Sep. 2013. 
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Figure 3. - Bayesian tree of the cytochrome c oxidase subunit I (COI) for 194 individuals of Gasterosteus spp. and other gasterosteids. The 
mean a posteriori values of the parameters are: TL = 38.329096; alpha = 0.076817; pinvar = 0.308267. As a reminder, G. aculeatus is full 
plated, whereas G. gymnurus has two to 10 lateral plates on both sides. 
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and reproduction data (Higuchi et al ., 2014). Besides, the 
plateless morph has been observed from Norway to Alaska 
via Iceland (Barrett, 2010), and is not characteristic to any 
European area. “Presence of plates only on pectoral region” 
is the sole character in the original description of G. gym¬ 
nurus. Moreover, as no type specimen is currently known 
(Kottelat, 1997), no taxonomical review on a reference indi¬ 
vidual is possible. 

The DNA taxonomy analysis with the COI marker 
does not differentiate G. aculeatus from G. gymnurus and 
G. islandicus. Geographically, a few divergences were 
observed on the mitochondrial cytochrome b and control 
region markers between European, American and Japanese 
populations, because of recent colonisations estimated at 
90,000 to 260,000 years from Alaska and British Columbia 
to Atlantic (Orti et al., 1994; Makinen and Merila, 2008). 
Several transatlantic exchanges have taken place during this 
period (Makinen and Merila, 2008). In addition, Makinen et 
al. (2006) demonstrate the absence of a relationship between 
phylogeny and habitat. While this lack of differentiation 
could be due to insufficient variability of the marker, it is 
here combined to the lack of morphological characters other 
than the highly plastic plate numbers. 

Mediterranean populations are distinct from West Euro¬ 
pean populations, however the dataset is not well representa¬ 
tive of the area (Makinen et al., 2006; Cano et al., 2008). 
Consequently, Bianco (2014) proposed that Mediterranean 
populations constitute another lineage. Our results also 
recover a Western Mediterranean subclade, but it does not 
include all Mediterranean specimens. Furthermore, the type 
locality of G. gymnurus is in England and North of France 
(Cuvier, 1829). 

In consequence, we consider Gasterosteus gymnurus 
Cuvier, 1829 as a junior synonym of G. aculeatus Linnaeus, 
1758. 

G. islandicus is endemic to Iceland, and is recognizable 
by its notch in pelvic girdle (Kottelat and Freyhof, 2007). The 
four specimens in our dataset are also included in the same 
cluster as G. aculeatus and G. gymnurus. As some popula¬ 
tions possibly belonging to this species were found in Nor¬ 
way and Finland (in sympatry with G. aculeatus ; Kottelat 
and Freyhof, 2007), the presence of the notch in pelvic girdle 
might also be a polymorphic character. More investigations 
are needed in order to clarify this species delineation. 
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